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Abstract. We summarize the dual gravity description for a thermal gauge theory, 
reviewing the key features of our holographic model of large N QCD [3l [4] and 
^ elaborating on some new results. The theory has matter in the fundamental 

JjJ ' representation and the gauge coupling runs logarithmically with energy scale at low 

, energies. At the highest energies the theory becomes approximately scale invariant, 

much like what we would expect for large N QCD although not with asymptotic 
^SJ . freedom. In this limit the theory has a gravity dual captured by an almost 

^ I classical supergravity description with a controlled quantum behavior, such that by 

renormalizing the supergravity action, we can compute the stress tensor of the dual 
, gauge theory. From the stress tensor we obtain the shear viscosity and the entropy 

' of the medium at a temperature T, and the violation of the bound for the viscosity 

00 ' to the entropy ratio is then investigated. By considering dynamics of open strings 

I in curved spacetime described by the supergravity limit, we compute the drag and 

diffusion coefficients for a heavy parton traversing the thermal medium. It is shown 
that both coefficients have a logarithmic dependence on momentum, consistent with 
pQCD expectations. Finally, we study the confinement/deconfinement mechanism 
' for quarks by analyzing open strings in the presence of the flavor seven branes. We 

flnd linear conflnement of quarks at low temperatures, while at high temperatures the 
quarkonium states melt, a behavior consistent with the existence of a deconflned phase. 
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1. Introduction 

A decade of running RHIC, tlie Relativistic Heavy Ion Collider at Brookliaven National 
Laboratory on Long Island, has produced a plethora of intriguing experimental results 
that have challenged our understanding of hot and dense strongly interacting systems. 
One of these remarkable findings is that the matter produced at RHIC appears to display 
striking flow characteristics that are consistent with expectations from ideal relativistic 
hydrodynamics In spite of the fact that the era of heavy ion physics at the LHC has 
barely begun, it appears that those exciting features also are present at higher energies 
[2]. The commonly accepted theoretical interpretation of this body of data is that 
the quark-gluon plasma (QGP) formed in the relativistic nuclear collisions performed 
at RHIC and the LHC is "strongly coupled": it displays aspects at odds with what 
perturbative QCD would predict. 

Strongly coupled quark-gluon plasma (sQGP) poses theoretically challenging yet 
experimentally accessible questions. The formation of QGP at RHIC is an example 
where theoretical descriptions are completely lacking at low energies because our 
perturbative techniques fail at strong couplings. However probing the non-perturbative 
regime of large N gauge theory through a gravity dual has led to some interesting 
results for the physics of the quark-gluon plasma. A popular approach so far has been 
the AdS/CFT correspondence [HIE], even though QCD is not a conformal field theory. 
However, for certain gauge theories with running couplings, there exist gravity duals. 
At zero temperatures these gravity duals are studied in [8] with, and in [7] without, 
fundamental flavors. On the other hand, at high temperatures, there are examples of 
gravity duals with or without fundamental flavors in the literature [TOl [TTl [121 [IS] ; all 
of which are quite successful in analyzing various aspects of strongly coupled nuclear 
matter in the IR. But each of the models have their limitations, most prominently due 
to the absence of UV completions, and a supergravity description for strongly coupled 
QCD incorporating various phases of quark matter is yet to be discovered. 

While the low energy limit of Klebanov-Strassler model with fundamental matter 
(which we call the Ouyang-Klebanov-Strassler or the OKS model) gives the most 
promising candidate to study large N strongly coupled QCD, the effective number of 
degrees of freedom grows indeflnitely in the UV. In the presence of fundamental flavors 
such a proliferation of degrees of freedom lead to Landau poles and UV divergences of 
the Wilson loops. However, since this cascading theory provides an ideal setup to study 
phase transitions as one can describe various degrees of freedom in a single framework, 
controlling the UV behavior then becomes essential to describe a consistent theory that 
could describe physics at all energies. In ^ such a UV completion was proposed and 
it is the goal of this paper to summarize and highlight some of the key features of this 
model and discuss its relevance to strongly coupled thermal QCD. 
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2. Gauge theory and the Dual Geometry 

Our goal is to construct the dual gravity of thermal field theory resembling features of 
strongly coupled QCD. A gauge theory which becomes almost conformal in the UV with 
no Landau poles or UV divergences of the Wilson loops, but has logarithmic running of 
coupling in the IR, can be thought of as a suitable candidate mimicking QCD. Following 
the developments in [7| and [8], in p] it was proposed that such a gauge theory indeed 
has a dual gravity. In fact the original proposal in [3] and |1] was more along the reverse 
direction, namely, the IR geometry of [8l was modified in steps in [1] to construct 
a dual gauge theory that may have the required properties to be in the same class as 
large QCD. On the gauge theory side we need IR confining and UV conformal. So 
the question is what should we do to the gravity dual at large r, i.e at large radial 
coordinate, to achieve the required properties? 

The answer, as discussed in details in [1] , is to first switch off the three- form fluxes 
at large r so as to keep only flve-form fluxes. This way AdS type solutions could be 
constructed. Secondly, to get rid of the Landau poles, we need to control the axio-dilaton 
behavior so that they decay rapidly at large r instead of having their usual growing 
logarithmic behavior. Both these were achieved by switching on anti D5-branes and 
anti seven-branes with electric and magnetic fluxes on them to kill off the unnecessary 
tachyons and thereby restoring stability. Additionally the anti seven-branes, alongwith 
the seven-branes, were embedded in a non-trivial way so as not to spoil the small r, or 
the IR, behavior of the theory [1]. The anti D5-branes are dissolved on the seven-branes, 
but they contribute to the three-form fluxes. Together they control both the three form 
fluxes as well as the axio-dilaton at large r. 

To be a bit more precise, the knowledge of the NS flux B2 and dilaton 
uniquely determines the running of the gauge couplings. On the other hand the dual 
geometry is described by the modifled conifold |1] and it is sourced by three form flux 
dB2 = H^s, Hrr, flve form flux and D5-D7 branes. In the UV, D5-branes source three 
form fluxes that exactly cancel the three form fluxes originally sourced by the wrapped 
D5 branes This way the effective three form fluxes vanish resulting in B2 = 0. Then 
we can write gi = g2 = Qym at large r and we get gyM ~ e*^/^ — thus the dilaton 
uniquely determines the gauge coupling. If r = Co + ie''^ is the axio-dilation fleld, 
then using F-theory [9] one can show that r behaves like an inverse power series for 
large r, i.e r ~ ^, z > [1], provided we have controlled the logarithmic behavior by 
anti-brane embeddings discussed above. In a more mathematical language of F-theory, 
r is determined by flnding the roots of the discriminant polynomial A(z) where z is the 
complexifled coordinate comprised of the two real coordinates of an internal P^. The 
roots give the precise location of the seven-branes and hence the dilaton is uniquely 
determined by the seven-brane embeddings. Using the inverse power series behavior of 
r, and identifying the energy scale with radial coordinate r ~ A, one gets for the gauge 
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coupling Qi = g2 = Qym the following behavior at the scale A: 

9ym ~ -J- (1) 

which indicates that the field theory reaches conformal fixed point at the highest 
energies. 

Motivated by the above construction in [3l H], the question now would be how to 
study the gauge theory, or the brane side of the story. The IR physics is captured by 
the near horizon geometry of D3 branes at the conifold point with additional wrapped 
D5-branes. Thus in this picture one needs to add the seven-branes as well as the anti- 
branes in addition to modifying the geometry. As we know, any brane configurations 
that are away from the original Klebanov-Strassler configuration, i.e not coinciding with 
the D3 and wrapped D5-branes, would survive on the dual gravity side as branes (and 
not convert to pure geometry and fiuxes). On the other hand, the modifications that 
we did on the gravity side to get to the required gauge theory may be too difficult to 
simulate using simple arrangements of branes. However we can still incorporate many 
of the informations on the gravity side to the brane side of the story, that could at least 
give us a hint of the required gauge theory. Thus a possible proposal would be that 
the following configuration of D3-D7-D5-branes (including the required anti-branes plus 
world-volume fiuxes) would give rise to our desired gauge theory, albeit in the large N 
limit with matter in fundamental representation: Place D3 branes at the tip of a 



N D3 Branes 




M Anti 5 Branes 



Figure 1. A very simplified brane configuration in ten dimensional geometry 
M.A X Cg, where A^4 is Minkowski space and Cg is a six dimensional cone. The actual 
configuration is non-trivial to depict pictorially and may not even exist as a simple 
arrangement of branes. The gravity dual, on the other hand, is more robust and can 
be easily described as shown in [4]. 

six dimensional conifold with radial coordinate r and base T^'^ = S*^ x S^. The world 
volume of the D3 branes extend in four Minkowski directions. Now place M number 
of D5 branes which wrap the S"^ on the conifold base and extend in four Minkowski 
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directions and place another M number of anti D5 branes away from but parallel to the 
D5 branes. The D5-D5 branes are separated along the radial r direction. Furthermore 
place Nf number of D7 branes extending in the radial r direction, wrapping the base 
T^'^ and filling up four Minkowski directions. The brane - anti brane configurations will 
result in tachyonic modes but introducing fluxes sourced by the branes can stabilize the 
system. Details of the embedding can be found in |3], [8] and one can, with the risk of 
oversimplification and possible pitfalls, depict the brane configuration as in Fig. 1. 

The excitations of these D3-D5-D7 branes (and the required anti-branes) are 
described by a gauge theory which in the UV has SU{N + M) x SU{N + M) color 
symmetry group and SU{Nf) x SU{Nf) flavor symmetry group [1] from the seven- 
branes. This is because addition of {p, q) branes at the junction, or more appropriately 
anti five-branes at the junction with gauge fluxes on its world-volume, tell us that the 
number of three-branes degrees of freedom are N + M, with the M factor coming from 
five-branes anti-five-branes pairs. Furthermore, the SU {N + M) x SU {N + M) gauge 
theory will tell us that the gravity dual is approximately AdS, but has RG flows because 
of the fundamental flavors. At the scale r = tq we expect one of the gauge group to be 
Higgsed, so that we are left with SU{N + M) x SU{N). Now both the gauge fields flow 
at different rates and give rise to the cascade that is slowed down by the Nf flavors. 
In the end, at far IR, we expect confinement at zero temperature. If gi,g2 are the two 
gauge couplings, they are related to the dilaton (p and the two form NS-NS flux B2 by 
following duality relations that capture the essence of gauge/gravity duality here: 
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A little bit of elaboration may elucidate the story further. As the anti five branes are 
placed at some radial location tq, their excitations give massive modes with energy scale 
Aq ~ tq. At low energies, A < Aq, these modes are not excited and hence can be ignored. 
This means at low energies, the effect of the anti five branes D5 will not be important 
and we effectively have Ouyang-Klebanov-Strassler geometry of [3j with a black-hole, 
which is region 1 of [4J. Thus in the IR we will have the SU{N + M) x SU{N) color 
symmetry group and SU{Nf) x SU{Nf) flavor symmetry group with gauge couplings 
5^1,5^2 running logarithmically as B2 ~ Mlog r. More precisely, using the exact form of 
fluxes and the dilaton at small r, i.e at low energies, we get the following beta functions: 
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At the IR, we see that the two gauge couplings run in opposite directions and SU (N+M) 
flows to strong coupling. Performing a Seiberg duality transformation, we identify the 
strongly coupled SU{N + M) with a weakly coupled SU{N - (M - Nf)) at the IR. We 
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see that not only the number of colors are reduced, but the difference of the size of the 
gauge group now decreases from M in to M — Nf. This difference will decrease by 
the increments of Nf until it is smaller than or equal to Nf. Then there are two possible 
end points: 

• If is still greater than zero then we will have an approximately conformal theory in 
the far IR, or 

• If decreases to zero but with finite M left over then we will have a SU{M) theory 
with Nf flavors that confines in the far IR (see [8] for more details). 

The latter theory, or more particularly the high temperature limit of the latter theory, 
is what we are interested in and henceforth we will only consider that theorjjjj. Once 
we know the dual background we can say that the Hilbert space of the gauge theory 
can be obtained from the Hilbert space of the string theory on this geometry. The dual 
geometry of above brane theory in the decoupling regime was proposed in One 
basically minimizes SUGRA action with the fluxes running logarithmically with r for 
small r and as inverse of r for large r. The resulting metric in the non extremal limit is 

yjh g[r) V/i 

= gapdx'^dx'^ + r'^Qmndx'^dx^ = Gf,„dx'^dx'' (4) 

where a,/3 = 0,1,2,3,4, m,n = 5,6,7,8,9 and /i, = 0,..,9. The internal five- 
dimensional metric is denoted by grnn whereas the other five- dimensional space is denoted 
by ga/B- We have also defined the following: 



^4 

h = h{r,'ip,(l)i,(l)2,0i,92) = -J 



J.4 



for large r (5) 
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i=l 

&jyfc(rft,V,^i,^2)log^(r) 



.j,k=0 



for small r 



for the warp factor h. Note that h is taken to be functions of all the internal as well 
as the radial coordinate r. This should be compared to earlier work of [7J where the 
warp factor is only a function of the radial coordinate. The rest of the variables are now 
defined as: 

H ^ H{r, 4^, 01, 02, e,, 92) = for large r (6) 

1=0 

= y dAr.,^^OMlog\r) ^^^^^ ^ 

j,k=0 

with (yf(r) = 1 — -j:!, where Vh is the black- hole horizon and Oj, bjk, Ci, djk are determined 
by the location of the D7 branes, the fluxes and the horizon. As discussed in [18], the 

I For a review of Seiberg dualities and cascading theories, see (TB]- For a review on brane constructions 
for cascading theories see [IT] . 
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non extremal solution has a regular Schwarzcliild horizon that covers the IR singularity 
of the extremal metric but the internal space gets modified. This modification of the 
internal space is described by the metric gmn which is not the metric of T^'^ but accounts 
for the deformation of the internal space. The forms of a^, bjk, Ci, djk and Qmn are given 
in |1] for the extremal case and the extremal solution of the IR geometry is discussed in 
[3]. The non extremal solution for the brane configuration in Fig 1 is under preparation 
[20] . However, the study of confinement mechanism and the computation of diffusion 
and drag coefficient will only require the knowledge of the warp factors at the location 
of the string which we will choose to be fixed in the internal directions. This means 
{ip, 01, 02, 01, G2) will be constant along the string and we will treat aj, hjk, Ci, djk as only 
functions of horizon radius Vh. 

To compute the stress tensor, we do need to know the exact form of Oj, bjk, Ci, djk 
as a function of internal coordinates. However, to obtain information about the four 
dimensional field theory, we will first integrate over the internal directions and after the 
integrations, we will be left with a five dimensional metric gai3 of the form given in (jlj) 
with warp factors h, H now only being functions of r. Then the coefficients Oj, bjk, Ci, djk 
of the five dimensional warp factors would only depend on r^. 

The question now is how does the warp factor h depend on the black-hole horizon 
r/i? Of course the Einstein equations along with the fiux equations will dictate this 
dependence- but one can restrict the form of the warp factor as a function of Vh by 
analyzing the dual gauge theory. In particular, the warp factor is related to the size of 
the dual gauge group [7] [3], 



at large r. Now from the T-dual description of the brane configuration in Fig 1 [19], 
one concludes that A^efr keeps on growing with energy scale A i.e. the radial coordinate 
r on the dual gravity side, which is depicted by Fig 2. From [7] we see that this means 



for large r. But a^'s are functions of Vh, thus ([H]) can be used to restrict the form a^s 
as functions of Vh- For example one way is that all ak{rh) < for all rh so that this will 
automatically satisfy (IHl). 

On the other hand for small r, one possibility is that the IR geometry is completely 
hidden behind the black-hole horizon, which will be the case for high temperatures [18j . 
Then we are left with the geometry at large r. However for low temperature and small 
r, we have 
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Again from the gauge theory interpretation that A^^efr grows with r, we get 
'jbjk log^r) kbjk log''~\r) 



E 
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(9) 



(10) 



As long as (E]) is satisfied at large r and ( TTOl) is satisfied at small r, we have a consistent 
supergravity description for thermal cascading gauge theory which fiows according to 
Fig 2. 



UV 



Scale A 



IR 



SU(N+M) XSU(N+M) 
Asymptotic AdS geometry 



Flow from UV to IR 



SU(N+M) X SU(N) 
g-log A 

Klebanov-Strassler type geometry 



Flow from UV to IR 



SU(M) 
g~Iog A 

Klebanov-Strassler type geometry 



Figure 2. The flow in gauge theory and its dual gravity description. A cascade of 
Seiberg dualities reduces the rank of the gauge group as energy scale is altered from 
UV to IR. 



3. The Stress Tensor and rj/s 

In order to compute expectation values of certain observables of the gauge theory, one 
first computes the partition function. As the Hilbert space of the quantum field theory 
arising from brane excitations in Fig 1 is contained in the Hilbert space of the geometry 
sourced by these branes, there should be a one-to-one correspondence between the 
partition functions of the gauge theory and the dual gravity. A precise mapping for 
the conformal case was proposed by Witten [6] who stated that the partition function 
of the strongly coupled quantum field theory should be identified with the partition 
function of the weakly coupled classical gravity. 

In particular if we are interested in computing the expectation value of an operator 
(O) with a source (po, one can make the following identification of the partition function 
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functional of the source 0o- 

>2gauge[0o] = (exp / 0oC) = 2gravity[0o] 

= exp(S'sUGRA[0o] + 5'gh[0o] + 5'counterterm[0o]) (H) 

where is a Minkowski manifold, S'gh is the Gibbons-Hawking boundary term , 0o 
should be understood as a fluctuation over a given configuration of field and 5'counterterm 
is the counter-term action added to renormalize the action. We will briefly describe how 
each term is obtained from supergravity action and their necessity. 

First observe that is a bulk field, which is in general a function of the coordinates 
of the geometry and 0o is the boundary value of 0. One integrates the classical ten 
dimensional gravity action over the compact manifold A^s to obtain an effective action 
^eff ^ dimensional manifold. The bulk five dimensional geometry has a boundary 
at r = oo and the boundary describes a four dimensional manifold with an induced 
Minkowski flat metric. Thus the boundary value of 

(f)o{t,x,y,z) = (f){r = oo,t,x,y,z) (12) 

is a field which lives in four dimensional flat space and plays the role of source in the 
four dimensional gauge theory. 

Now to obtain the gravity action as a functional of (po, one integrates the five 
dimensional effective action Sf^ over the radial coordinate r and we are left with gravity 
action S'[0o] as a function of 0o- It turns out that in general for the gravity action S[(j)] 
to be stationary under perturbation — > + 50, one needs to add surface terms S'gh[0] 
which are known as the Gibbons- Hawking terms [22]. On the other hand, after the radial 
integral is done, for some sources 0, S'[0o] becomes infinite and one needs to regularize 
the action to obtain finite expectation values for operators. The regularization requires 
addition of extra terms which cancel the infinities that appear in the gravity action and 
are denoted by 5'counterterm- Keeping all this in mind, one takes the functional derivative 
of (fTTl) with respect to 0o to obtain 

_ ^^gauge _ ^^gravity 
500 500 

~ exp (^total 00 ) 7- (13) 

500 <f>o=0 

where we depict S'totai[0o] as: S'totai[0o] = 5'[0o] + 5'gh[0o] + 5'countcrtcrm[0o]- Note that, 
from our construction, this is also the renormalized action. 

Observe that in the usual AdS/CFT case we consider the action at the boundary to 
map it directly to the dual gauge theory side. For general gauge/gravity dualities, there 
are many possibilities of defining different gauge theories at the boundary depending 
on how we cut-off the geometry and add UV caps. The details on adding geometries at 
large r to existing IR geometries for small r can be found in |31 H]. 

One of the important thermodynamical quantity to extract from our background 
would be the stress tensor. As the stress tensor couples to the metric, its expectation 
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value is obtained from the renormalized action S'totai via (T*-^(Ac)) = ^'^^'l?'^''^^ 

where we have chosen 0o = hj as the four dimensional effective metric induced by 
the ten dimensional metric of the form (jll) (see [3] for details on the procedure). At 
every fixed r the five dimensional metric induces a four dimensional metric and with 
appropriate rescaling of the time coordinate, we can view this effective four dimensional 
metric as a fiat Minkowski metric (with possible perturbations). If Oij denote the 
perturbations on the background by quark strings that stretch from the D7 branes to 
the horizon of the black hole, then writing the supergravity action up to quadratic order 
in Oij] using integration by parts together with appropriate Gibbons-Hawking terms; 
and holographically renormalizing the subsequent action, we obtain: 



rpmm f ^ ^ \ ^ / (' W™" _L Zl7'"-'"'\ r.(4)[/3] /I (' JV'™-" i T<r^'"^\ o 

-'mcdium+quark — y (27r)'^ ^ " \ ^ ~^ '^\a\ J^'nn ^\'^\a\ "r-'^|a|J'^ 



(4)[/3] 
nn 



+ (X|7 + irfjp).(f^[''l + J"5„^e-^--^- + OiTe-^-^) (14) 

where (^^'(]) , A/'uv) together will specify the full boundary theory for a specific UV 
complete theorj|§|. The quantities jj^^^j^^x^^ etc. with m,n = 0,..,3 are given 
in [H H] and are completely independent of the radial coordinate r and depend on the 
3+1 dimensional spacetime coordinates. Note also that the stress tensor depends only 
on the temperature T and is independent of any cut-offs. The dependence on the UV 
degrees of freedom is exponentially suppressed, and in the limit Af^y = e~",n >> 1 we 
reproduce the result for the parent cascading theory. 

This exponential suppression comes from the form of A^^efr = r'^h where h is the 
warp factor. For geometries such that N^s^r) ~ A^(l + log r) at large r, one easily gets 
that ~ gi^cffC'") = giA/'uv -^Iiqyq have taken r — )■ oo limit and defined iVeff(oo) = Af^v 
Sending r — )■ oo also means Afuv — ^ oo, thus the final result for the stress tensor will 
be same for various UV completions. However, the point we are trying to make is that 
different UV completions reach the same IR limit. 

For Asymptotic AdS spaces, we have the warp factors behaving like ([5]) and 
Ncs{oo) = Afuv = N{1 + ao). Then for the stress tensor we get 



(4)[/3] 
nn 



rpmm _ f ^ 1 \ ^ j (TTmn . TTnm\ {4:)[I3] AfT^mn . K^mx 

-'medium+quark — y ^27r)^ \ l"!'^ "" ^\'^\a\ "'"-"■|a|J'' 

+ (i^i^p + irrjp).(f2['^] + (15) 

where an°^ is uniquely determined by Oq which appears in (|9]) and thus determines Muv 
Here we have already set terms to zero, so they don't appear in ( IT5|) . and dependence 

§ It may appear from ()14p above that the energy-momentum tensor is diagonal. This is however not 
so because we are using the variables defined in cq. (3.108) of [3], to write However if 

we replace Lmm by Z^i/, defined in eq. (3.107) of j^, the energy- momentum tensor will not appear 
diagonal. We simply found it convenient in [3l to express the energy-momentum tensor using the Lmm 
variables. These subtleties are described in full details in [3]. 
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on Muv appears implicitly through ai"'' . More precisely, for a generic asymptotic AdS UV 
completion, there might not be any analytic inverse functions T such that r ~ J^{Muv)- 
Hence it might not be possible to see how ~ J^'^Muv) reaches zero. But the final 
result is independent of how this limit is reached and is only sensitive to ai"^ , thus (|T5l) 
uniquely determines the field theory stress tensor. 

With the general formulation of stress tensor (fT^.( |T5|) . which is similar to the 
AdS/CFT results [21] in the limiting case M = iVj = 0, we can compute the wake 
a moving quark creates in a plasma with Oij being the metric perturbation due to 
string (see |23] for an equivalent AdS calculation). Furthermore we can compute the 
shear viscosity r] from the Kubo formula with the propagator obtained from the dual 
partition function ( |TT|) and the entropy density s obtained by computing energy and 
pressure density using the stress tensors (JH]), (fT5|) . Alternatively entropy s can be 
identified with the black hole entropy using Wald's formula and the result for the ratio 
77/s is given by: 

Tj 1 

s 



47r + i log2(l-r4e-4A/'-) 



(16) 



C3K 



3L2 (1 - T^e-^- 



n3/2 



5,(47r2-log2C,)+47rA log C„ 



47r2 - log^ Co + 167r2 log^ Co 



< 



where {Ao, Bo, Co, ak) are constants that depend on the temperature T and e~-^"^; and 
C3 is the coefficient of the Riemann square term coming from the back reactions of the 
D7 branes in the background. These have been explicitly worked out in [3]. The C3 
dependence of the rj/s ratio first appeared in |2l] and the above expression shows that 
for C3 > 0, we have violation of the lower bound for rj/s. For C3 greater than some critical 
value, studies have shown that the theory becomes acausal [25]- [32]. This restricts the 
form of the terms that appears from D7 brane back reaction as we require the boundary 
gauge theory to be causal. Observe that for asymptotic AdS space, using (fTSll for stress 
tensor, rj/s is given by f|T6|) with Af^v = 00 and then the only violation comes from the 
Riemann square term. 
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4. Momentum broadening and drag 

Consider a parton moving through a plasma in four dimensional Minkowski spacetime 
with the following world line 

x{t) = vt + 6x(t) 

z{t) = y{t) ^ 5y{t) (17) 

For a fast moving parton, we can always choose coordinates such that f|T7j) is the world 
line. Now if the plasma has matter in fundamental representation, has logarithmic 
running of coupling in the IR that becomes asymptotically conformal, we can treat the 
ten dimensional geometry with metric (jl]) to be the dual gravity of this gauge theory 
which lives in four dimensional fiat spacetime. 
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To obtain the momentum broadening of the parton we shall use the Wigner 
distribution function / as defined in QCD kinetic theory [31] 

/(X,rx) ^ (/cc(X,rx)) = Tr [pQt(X^)Uab(X^, Xi)Qa(Xi)] 

/(X,ri) ^ {fcc{X,rL)) = Tr [pQKX^ )Uab(X^, )Q,(X^)] (18) 

where = X - r±/2,X^ = X + r±/2,X^ = X - rL/2, X^ = X + rL/2, X = (t, x) 
is the world line of the parton field Qa without any fluctuation, Uab is the link and p 
is the density matrix. The index ab refers to color and for our choice of the world line 
for the parton we have = 25y'^,r\ = dx^. Now Fourier transforming the distribution 
functions, we can get the average transverse momentum to be 

{pD = jd'^j ^^klfiXM) (19) 
whereas the average longitudinal momentum fluctuation is given by 

(pi) = fd'x [^klf{X,k,) (20) 



. . (27r)2 

Now if we assume that the initial transverse momentum and longitudinal momentum 
fluctuation distributions are narrow, then we have 

{pD =2^Tr 

{pD =^lT (21) 

with T some large time interval, and kt, are the diffusion coefficients. Generalizing 
the arguments in [35] to include longitudinal momentum fluctuations, one can write the 
diffusion coefficients solely in terms of functional derivative of Wilson loops. The final 
result is 

kt = hm ^ / dt e^^* [iG\^{t, 0) + iGl^{t, 0) + iG\^{t, 0) + iGl^{t, 0)] 
w— ^0 4 J 

kl = lim ^ / dt e*^* [iGu{t, 0) + iGl^{t, 0) + iGl^{t, 0) + iGl^{t, 0)] (22) 
i^— >o 4 J 

where the Green's functions are 



trp»Wc[0,01\ ^Ci(t)i5Ci(t') 



trpOWc[0,0]\ ^ 5C2(t)<5C2(t' 

G^ it t') - 1 /t„o ^^Wc[^Ci,^C2 
'^^^^*'*^-trpOWc[0,0]\''^ 5Ci(t)5C2(tO 

with ( = y,x and t, t' are the real part of complex time tc, t'(j on contour C (fig 3). Here 
we have introduced type '1' and '2' fields i = 1, 2) of thermal field theory evaluated 
along the upper and lower horizontal line of the contour G. We denote by VFcli^Ci) ^C2] 
the Wilson loop with deformation 5Ci and 5(^2 on the upper and lower lines of G . 
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Xi=(t„,x„ ,y(t„)) 




X=(t,x +vAt,y(t)) 



X,. =(to-is,}(o-is,y(to-is)) 



Figure 3. The contour C with X denoting the coordinates of a point on the contour. 
The upper line with real time coordinate is the worldline of the heavy parton while 
the lower line has complex time coordinate. Fields evaluated on the upper line are the 
type '1' fields while those evaluated at lower line are type '2' fields of thermal field 
theory. Here s — >■ is real. 



The Wilson loop at strong coupling can be computed using holography. This means 
we identify 

{tip^Wc) =e^^NG (24) 

where p° is the density matrix [35], S'ng being the Nambu-Goto action, with the 
boundary of the string worldsheet being the curve C. That is the string worldsheet 
ends on the world line f|T7j) of the heavy parton in four dimensional flat spacetime. If 
: (cr, r) — )■ {t, x,y,z,( = ^/r, ip, 0i, 02, 6*1, O2) is a mapping from string worldsheet to 
ten-dimensional geometry, we restrict to the case when 

X° = t, X^ = r, X^ = x{T,a), = = 6y{T,a) 

X" = {k = 5,Q,7), {X\X') = (61,62) =71 (25) 

With the choice of parameterizatioijjj] t = t,a = C, also known as the static gauge, we 
have 

x{t, r) = vt + x{r) + 5x{t, r) 

z{t,r) = y{t,r) = 6y{t,r) (26) 

where x is the unperturbed solution for the mapping. Minimizing the Nambu-Goto 
action gives the fields x(t,r),y(t,r) and z(t,r) and finally integrating S'ng over radial 
direction r, one gets the boundary action. We can obtain the Green's function from the 
boundary action and the result for the diffusion coefficients are [39] 

KT = 7r^jgsN,sT^{l + B) 

= 7'/'vrv/^JVe^r3(l + B) (27) 

where 7 = l/Vl — "W^, i3 is a constant of 0{gsM'^ /N, g'^M'^Nj /N) which depends on the 
velocity v, and is defined in the following way: 

f log(7) - b log(r;,)log(7) + |log2(7) 



B 



1 - a hg{rh) + b log^(r/,) 
3gsM^ 3g!M^Nf 



(28) 



II We denote the axio-dilaton and the world-sheet time using the same symbol r (as is the usual lore 
in string theory). It should be clear from the context which one is meant. 
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Here N^g is the number of effective degrees of freedom for the boundary gauge theory 



327r2A^ 47riV V 47riV 167r2A^ / V 



(29) 



where is the radial distance where the warp factor changes from a logarithm to a 
power series and tq is associated with the maximum depth of the seven branes in the 
radial direction. The value of and tq will of course be determined by the embedding 
equations for the branes which give their precise location in ten dimensional space. Thus 
for different configuration of branes in ten dimensions, we will end up with different 
effective degrees of freedom for the gauge theory in four-dimensional Minkowski space. 
Note for duality to hold, must be quite large, making N^^i quite large. In deriving 
(12^ above, we utilized the logarithmic behavior of the warp factors h, H for r and 
the precise form for h given by. 



(30) 



2tiN 

For A^j = M = 0, we get back the value of kt,hl as computed in [3S][5nj. However 
for a non conformal field theory with fundamental matter - which is more relevant for 
QCD - M 7^ 0, Nf 7^ 0, and our analysis shows that the diffusion coefficients kt and kl 
depends non trivially on velocity v of the parton. In fact kt, kl varies as log(l — v^) 
and increases with velocity. This dependence on velocity is consistent with perturbative 
calculation of diffusion for heavy partons due to collisions with quarks and gluons using 
Langevin dynamics [38j. Details of our calculation can be found in [39] . 

Now we compute the drag experienced by a parton as it ploughs through a thermal 
plasma at constant speed using dual gravity description where the parton is represented 
by a string moving in ten dimensional geometry with metric (jlj). We consider the map 
( 125|) . (126|) in the limit 6x = 6y = 0. From the Nambu-Goto action one derives the 
following canonical momenta for the string: 



2/ 



m . (31) 



where Tq is the string tension l/27ra', Gij is the metric given by (jl]), /C is the determinant 
of the world sheet without any fluctuations, prime denotes derivative with respect to r 
while dot denotes derivative with respect to t. Now the rate at which momentum is lost 
to the black hole is given by the momentum density at horizon 

Tll{r = rn) = -T,Cv (32) 

while the force the parton experiences due to friction with the plasma is ^ = —up with 
p = mv/ \/l — v^. To keep the parton moving at constant velocity, an external field Si 
does work and the equivalent energy is dumped into the medium [1Q][1T]. Thus the rate 
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at which a quark dumps energy and momentum into the thermal medium is precisely 
the rate at which the string loses energy and momentum to the black hole. Thus up to 
0{gsNf,gM) we have ^ = -n^(r = rn) and 



27ia'm 
1 



(33) 



h 



2na'mL'^ 



2-kN 



log 



(1-1)2)1/4 



1 + 



2-n 



log 



+ 



rj^2 



(34) 

log(7), similar to B 



2'Ka'm 

which defines C and it depends on velocity v. More precisely, C 
given in fl28l) . 

In determining we used the logarithmic behavior for the warp factor h near 
horizon r = rh given by (l30l) . The above result for drag should now be compared with 
the AdS result [10] [H]. In the AdS limit, M = A^/ = and we get back the AdS 
drag. However, unlike the AdS case where drag is independent of velocity, the result in 
( l33l) is velocity dependent. In particular it decreases with increasing velocity- a result 
consistent with the perturbative analysis of [38] . 



5. Confinement /deconfinement from holography 

The linear confinement of quarks at large separations and low temperatures is a strong 
coupling phenomenon and one uses lattice QCD to compute the free energy of the bound 
state of quarks [121 |l3l HI] . On the other hand at high enough temperatures, if the gauge 
coupling is weak, one can compute the free energy using perturbative QCD and obtains 
Coulombic interactions between the quarks. Whether using perturbative QCD at weak 
coupling or lattice QCD and effective field theory techniques at strong couplings, one 
finds that at high temperatures, the Coulomb potential is Debye screened [15] [i3] . 

The study of heavy quark potential gathers special attention as it is linked to a 
possible signal from Quark Gluon Plasma. In particular, heavy quarks are formed during 
early stages of a heavy ion collision as only then there exists enough energy for their 
formation. On the other hand right after the collision QGP is formed as only then the 
temperature is high enough for quarks to be deconfined. As temperature goes down, 
heavy quarks form bound states and as these states are very massive, they are formed at 
temperatures higher than deconfinement value. This means heavy quark bound states 
like J/iIj can coexist with QGP and can act as probes to the medium. 

In particular one can study the J / ip bound state formed in proton-proton collisions 
and compare with heavy ion collision where a medium is formed. The medium will screen 
the cc interaction making them less bound and eventually resulting in a suppression of 
J/ip production. This phenomenon is known as the J /if) suppression and considered as 
a signal of QGP formation [16]. 



16 



In order to quantify quarkonium suppression and analyze features of plasma that 
causes the suppression, one has to compute free energy of J / ip. At large couplings, lattice 
QCD calculations are most reliable and there has been extensive studies quantifying the 
potential for heavy quarks |17]-[56] . The effective potential between the quark anti- 
quark pairs separated by a distance d at temperature T can then be expressed succinctly 
in terms of the free energy F{d,T), which generically takes the following form: 

F{d,T) = adfMT)-^fcid,T) (35) 

d 

where a is the string tension, a is the gauge coupling and fc and fg are the screening 
function^ (see for example [lT]-[53] and references therein). At zero temperature free 
energy is the potential energy of the pair. On the other hand free energy or potential 
energy of quarkonium is related to the Wilson loop which we will elaborate here. 

Consider the Wilson loop of a rectangular path C with space like width d and time 
like length T. Note that we are working with Euclidean metric with time t = it and 
the time like paths can be thought of as world lines of pair of quarks QQ separated 
by a spatial distance d. Studying the expectation value of the Wilson loop in the limit 
T — )■ oo, one can show that it behaves as 

(WiC)) ~ expi-TEQQ) (36) 
where Eqq is the energy of the QQ pair which we can identify with their potential energy 
YgQ^d) as the quarks are static. At this point we can use the principle of holography [5] 
^33j and identify the expectation value of the Wilson loop with the exponential of the 
renormalized Nambu-Goto action, 

{W{C)) ~ eM-S'^o) (37) 
with the understanding that C is now the boundary of string world sheet. Note that 
we are computing Wilson loop of gauge theory living on fiat four dimensional space- 
time Whereas the string world sheet is embedded in curved five-dimensional 
manifold with coordinates x^'^''^'^ and r. We will identify the five- dimensional manifold 
with Region 3 of [1]. For the correspondence in fl37|) to be valid, we need the t'Hooft 
coupling which is the gauge coupling in the theory to be large. On the other hand as 
discussed before, it is in this regime of strong coupling that linear confinement is realized 
in gauge theories. Thus using gauge/gravity duality is most appropriate in computing 
quarkonium potential. Comparing (l36l) and (1371) we can read off the potential 

Cren 

V,, = ^..J-f (38) 

Thus knowing the renormalized string world sheet action, we can compute Vqq for a 
strongly coupled gauge theory. 

For non-zero temperature, the free energy is related to the Wilson lines W (±|) 

via: 

F(rf,T)" 



exp 



T 



(l^t(+|))(iy(-|)) ^^^^ 



^ We expect the screening functions /s,/c to equal identity when the temperature goes to zero. This 
gives the zero temperature CorneU potential. 
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In terms of Wilson loop, the free energy fl35p is now related to the renormalized Nambu- 
Goto action for the string on a background with a black-holfl One may also note that 
the theory we get is a four- dimensional theory compactified on a circle in Euclideanised 
version and not a three-dimensional theory. 

We will restrict to cascading gauge theories where the effective number of colors 
grows as the scale grows. This property of a gauge theory is most relevant for physical 
theories as new degrees of freedom emerge at UV and effective degrees of freedom shrink 
in IR to form condensates at low energy [16] . The number of colors at any scale u = 1/r 
is given by (|9]) and for the analysis given here, it is simpler to define 



instead of N^i^iu). In terms of 'H{u), the condition that Nqs^u) is a decreasing function 
of M = 1/r becomes 

n'{u) > (41) 

Combining fHOj) and fHTj) yields the following condition 

niu) > ^ (42) 
From now on, the value of L is set to 1 for the rest of this subsection, so that ?/(«) > 1. 

5.1. Zero temperature 

Let Mmax be the maximum value of u for the string between the quark and the anti- 
quark. Then the relationship between Wmax and the distance between the quark and the 
anti-quark is given by [1] 



d{u^a^) = 2Mjnaxrt(Umax) / dv —p^^ ^ 



'H(Mmaxt^) 



while the renormalized Nambu-Goto action can be written as [1] 

'T Umax t ~ -L Jq U 



-1/2 

(43) 



^ dv 



HiUr 



1 -1/2 

+ 0{eo) \ (44) 



where eo gives the depth of the seven brane on which the string is attached, Qm, Qm 
depends on the metric and is explicitely given in [4j for our supergravity background. 

+ There is a big literature on the subject where quark anti-quark potential has been computed 
using various different approaches like pNRQCD [M] — [SS], hard wall AdS/CFT [57] — [SI] and other 
techniques [62j— [66j . Its reassuring to note that the results that we get using our newly constructed 
background matches very well with the results presented in the above references. This tells us that 
despite the large N nature there is an underlying universal behavior of the confining potential. 
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Observe that for Wmax small, d is small and ignoring higher order terms in Mmax, we 
get from fHSl) and 



n 0*^^ 

= -— — + (0.174^2 + 0.095A)rf + O(d^) (45) 
d 

where A2, G2 depend on the geometry and their values for our metric dl]),© will depend 
on the coefficients a^, Cj which are computed in [4J. However, the potential is dominated 
by the inverse d behavior, i.e the expected Coulombic behavior and the coefficient of the 
Coulomb term which is a dimensionless number is independent of the warp factor and 
therefore should be universal. This result, in appropriate units, is of the same order of 
magnitude as the real Coulombic term obtained by comparing with Charmonium spectra 
as first modeled by [67] and subsequently by several authors [17]-[53].[M]-[56].[58]-[6T]. 
This prediction, along with the overall minus sign, should be regarded as a success of 
our model (see also [68] where somewhat similar results have been derived in a string 
theory inspired model). The second term on the other hand is model dependent, and 
vanishes in the pure AdS background. 

Now we will analyze ( H3|) and (jH]) for large Umax- First note that for (H3|) to 
represents the physical distance between a quark and an anti-quark in vacuum, the 
integral must be real. This is guaranteed if for all < f < 1: 

^(^^Kax)-^^^f:^^^) <1 (46) 

For AdS space, ( l46l) is automatic, as "H = 1 and then d is proportional to "Umax which 
results in only Coulomb potential for all values of d. But for a generic warp factor, (146|1 
gives rise to a finite upper bound for Wmax and we will analyze the behavior of d and the 
action S^q near the upper bound. 

To show confinement at large distances the potential between the quark and the 
anti-quark must be long ranged. That is, d{umax) must range from to 00 as Wmax 
varies from to it's upper bound, say Wmax = a;max- Since T-Li^u) > 1 for our non-AdS 
geometry, the only way to satisfy these conditions is via sufficiently fast vanishing of the 
square- root in (H3l) as f — )■ 1 at Umax = a^max- Studying the behavior of the square-root 
[1], one concludes that Xmax must be the smallest positive solution of 

xn\x) - 2n{x) = (47) 

For our dual geometry with metric (!4]),([5]), indeed (H71) has a real positive solution Xmax- 
As Mmax approachcs Xmax, both d and S^q become infinite, growing linearly with Wmax 
and we get [1] 

v., ^ (IH) " (48) 
i.e linear confinement of quarks separated by large distance at zero temperature. 
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5.2. Finite temperature 

At finite temperature, the relation between Wmax and the distance between the quark 
and the anti-quark is obtained by replacing 'H(-u) with a/1 — u'^/uf^'H{u) in 



/ v"^ /v. 

driu^ax) = 2Mmax\/ 1 " U^^J ulH{u„,a^) / dv — — 

Jen [l — V Mjnax/' 



4 (l-<ax/0 ( 

(1 - v^uijut) \n{u. 



eo V- - -max/0(^(Wmax^^))^ 
1/2 



(49) 



The explicit factor of Wmax makes ^^(^max) vanish at Wmax = as in the T = case. As 
"Umax approaches Uh, the integral near v = 1 behaves like 

driu^..) - f dv -<axK ^ 

Jo [I - V)[l - VUma.^/ Uh) 

which indicates that (ir(-Umax) goes to as Wmax approaches Uh- Hence, at both Wmax = 
and Mmax = dTiu^gji) vanishes. Since dxiu^g,^ is positive in general, there has to be 
a maximum between Umax = and Wmax = Uh- Whether the maximum value of (ir('Umax) 
is infinite as in the T = case depends on the temperature (equivalently, u^^) as we 
now show. 

The fact that the physical distance needs to be real yields the following condition. 
For all < < 1, 



W^t(.Iw) - ( ^P^^ < 1 (51) 



Similarly to the T = case, c^t (^max) can have an infinite range if Umax = 2/max is the 
smallest positive solution of the following equation 

yU\y) - 2U{y) = {yjuhf yU\y) (52) 

Note that the left hand side is the same as the zero temperature condition, ( 147|) . 
The right hand side is the temperature {uh) dependent part. Using the facts that: 

j—l oo 

JI (fc - 1/2) = -{23 - 3)!!/2^ J^x^' (2j - 3)!!/2^j! = 1 - Vl^^ 

k=Q j=l 



it is clear that y = Uh cannot be a solution of fl52|) because at y = Uh, the equation 
reduces to l-iiuh) = which is inconsistent with the fact that T-L{y) > 1. 

Recall that we are considering gauge theories for which ^{y) > 1 and 'H'{y) > 0, 
and thus the right hand side of (!52ll is always non-negative. As y increases from towards 
Xmax; the Icft hand side of 0521) increases from —2 while the right hand side increases from 
0. Assuming that the warp factors do not get corrections from black hole horizon, i.e. 
'H{u,Uh = 0) = l-iiu^Uh) for all Uh, we see that left hand side reaches when y = Xmax 
which is the point where the distance d{u^ax) at T = becomes infinite. At this point 
the right hand side of (j52l) is positive and has the value (a;max/^/i)'^a^max'H'(a;max)- Hence 
the solution of (l52l) . if it exists, must be larger than Xmax- 
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Even if we consider to be a function of horizon Uh, we expect solution of fl52p 
y = Umax > a:^max- This is because black holes will attract the U shaped string toward the 
horizon and a string attached to quarks separated by a distance d will stretch deeper 
into the bulk geometry in the presence of a black hole. A warp factor h{u, Uh) that 
explicitely depends on horizon and thus 'H{u,Uh) explicitely depending on Uh will not 
alter this attraction of the black hole and we will have ?/max > a;max- 

Consider first low enough temperatures so that Uh ^ a^max- For these low 
temperatures, (152|) will have a solution, as the right hand side will be still small around 
y = a^max- This then implies that the linear potential at low temperature will have an 
infinite range if the zero temperature potential has an infinite range and we get 



Now we show that the infinite range potential cannot be maintained at all 
temperatures. We can have a black hole such that Uh = Xmax- When the left hand 
side vanishes at y = Xmax, the right hand side is Xma.x'H' {x^aax) = 2'H(a;max) which is 
positive and finite. For y > Xraa.x, the left hand side (yH^y) — 2'H{y)) may become 
positive, but it is always smaller than yT-L'^y) since T-Li^y) is always positive. But for the 
same y, the right hand side {{y / UhYyT-i' {y)) is always positive and necessarily larger 
than yl-Ciy) since (y/uh) > 1. Hence, fl52l) cannot have a real and positive solution 
when Uh = Xmax- Therefore between Uh = oo and Uh = Xmax, there must be a point 
Uh = Uc when ( f5^) cease to have a positive solution. Thus for some Uh < Uc and hence 
for temperature T > Tc, ( 152|) has no solution and c^t (wmax) will not diverge for any Wmax 
within (0,M/i). Furthermore, since the expression vanishes at both ends, there must 
be a maximum (iT(Mmax) at a non-zero Wmax- When the distance between the quark 
and the anti-quark is greater than this maximum distance, there can no longer be a 
string connecting the quark and the anti-quark. Hence for T > T^. the string breaks 
and we have two parallel straight strings describing the quarks. Of course there will be 
interactions between the two straight strings but these integrations will be a' suppressed 
and thus can be ignored in the a' — )■ limit which is the decoupling limit where duality 
is precise. Hence we can treat these quarks described by straight strings as 'free' quarks. 

To summarize, we have just shown that if we start with a dual geometry that allows 
infinite range linear potential at zero temperature, there exists some critical temperature 
above which the string connecting the quarks break. This shows that at high enough 
temperatures quarkonium state melts and gives rise to 'free quarks'. 

Thus our dual gravity analysis indicates that if we have long range interactions of 
QQ pairs at zero temperature where the potential energy grows linearly with distance, 
this interaction cannot be sustained at higher temperatures. This is the key point of our 
analysis, as this shows that quarks that are linearly confined at low temperatures, must 
be deconfined at higher temperatures and thus we have a confinement/deconfinement 
phase transition. 

In Fig 4 we plot of interquark separation (i as a function of Wmax while in Fig 5 
we plot free energy denoted by Vqq as a function of separation d. We have identified 




(54) 
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the renormalized Nambu-Goto action in the presence of a black hole with the free 
energy of QQ pair at finite temperature i.e. F{d,T) ~ S^q/T where T = j dr with r 
being Euclidean time. This way of obtaining free energy is indeed consistent with the 
identification of partition functions of gauge theory with that of gravity. For the plots 
in Fig 4 and Fig 5, we have made a particular choice for the warp factors 

h = L\%xp{-au^), H = expipJ) (55) 

where {a, (3) = (3,3) and (a,/3) = (0.1,0.05). This choice is consistent with our ansatz 
(jlj), ([5]) and by arranging the D5 — D5 — D7 in similar fashion as in Fig 1, the warp 
factor in (!55l) can be obtained. 
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Figure 4. Inter quark distance as a function of Mmax for various temperatures and 
warp factor with (a, a, ^, fi) = (0.1, 3, 0.05, 3) in the warp factor equation. 



Observe that for a wide range of temperatures < 1/uh = T < l/uc = Tc, 
the potential and thus the free energy hardly changes. But near a narrow range of 
temperatures Tc — e<T<Tc + e (where e ~ 0.05), free energy changes significantly 
.Our numerical analysis suggest the presence of a deconfinement transition, where for 
a narrow range of temperatures 0.28 < < 0.39 the free energy of QQ pair shows a 
sharp decline. Interestingly, changing the powers of u in the exponential changes the 
range of Tc only by a small amount. So effectively Tc lies in the range 0.2 < Tc < 0.4. 
Putting back units, and defining the boundary temperatur^ T as T = — ^ ^"''^ , our 

* See sec. (3.1) of j3j for details. 
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Figure 5. Heavy quark potential Vqq as a function of quark separation d with cubic 
warp factor, or equivalently, (a, 5,/3,/3) = (0.1,3,0.05,3) in the warp factor equation 
for various temperatures. 



analysis reveal: 

0.91 - 1.06 

-JT <Tc< ^ (56) 

which is the range of the melting temperatures in these class of theories for heavy 
quarkonium states. Since the temperatures at both ends do not differ very much, this 
tells us that the melting temperature is inversely related to the asymptotic AdS radius 
in large thermal QCD. 



6. Conclusions 



In this note we have summarized our proposal for the gravity dual of a non conformal 
finite temperature field theory with matter in fundamental representation. The gauge 
couplings run logarithmically in the IR while in the UV they become almost constant and 
the theory approaches conformal fixed point. To our knowledge, the brane construction 
and the dual geometry presented here is the first attempt to UV complete a Klebanov- 
Strassler type gauge theory with an asymptotically conformal field theory. 

Although our construction is rather technical with the gauge group being of the form 
SU (N + M) X SU (N) in the IR, one can perform a cascade of Seiberg dualities to obtain 
the group SU{M) and identify this with strongly coupled QCD. One may even interpret 
that the gauge groups depicted in Fig 2 contain strongly coupled large N QCD. This is 
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indeed consistent as the coupling g^+M of SU{N+M) factor in SU{N+M)xSU{N+M) 
or in SU{N + M) x SU{N), always decreases as scale is increased. To see this, observe 
that in the IR, the coupling Qn+m runs logarithmically with scale and gets stronger as 
scale is decreased. At the UV, we can arrange the sources in the dual geometry such 
that Qn+m decreases as the scale grows and runs as Qm+n ~ Sfc/A'^. By demandinj^ 

> (57) 
we see that Qn+m indeed decreases as scale is increased. Thus from UV to IR coupling 
always increases, just like QCD. 

Of course the 't Hooft coupling Aiv+Af = (A^ + M)gN+M for the group SU{N + M) 
is still large in the limit N ^ oo, even if Qn+m has decreased to very small value. This 
allows us to use the classical dual gravity description for the gauge theory which has 
a coupling that shrinks in the UV, mimicking QCD. For even higher energies, 't Hooft 
coupling will eventually become too small for finite N + M and supergravity description 
will no longer hold. But we can use perturbative methods to analyze the field theory 
in the highest energies. For finite N + M and very high energies, the gauge coupling 
may even vanish giving rise to asymptotically free theory. These arguments lead us to 
conclude that, in principle, the brane configuration we proposed can incorporate QCD 
and in the large N limit, the gravity dual we constructed can capture features of QCD. 

Using the dual geometry, we have studied the dynamics of 'quarks' in the gauge 
theory and computed shear viscosity rj and its ratio to entropy rj/s. The key to most 
of our analysis was the calculation of the stress tensor of gauge theory and we showed 
how different UV completions contribute to its expectation value. Using the correlation 
function for the stress tensor, we computed the shear viscosity of the medium while the 
computation of pressure and energy density allowed us to calculate the entropy of the 
system. Using a similar procedure to calculate correlators of stress energy tensors, with 
introducing diagonal perturbations in the background metric, we can easily evaluate 
the bulk viscosity ( of the non-conformal fluid. One can consider vector and tensor 
fields of higher rank to couple to the graviton perturbations in the five dimensional 
effective theory and study how this coupling effects the bulk viscosity. The calculation 
is underway and we hope to report on it in the near future. 

All the calculations we performed regarding properties of the plasma did not account 
the effect of expansion of the medium which is crucial in analyzing fluid dynamics. One 
possible improvement would be to construct a time dependent dual gravity which can 
describe the expansion of the QGP formed in heavy ion collisions. A first attempt 
would be to consider collisions of open strings ending on D7 branes and then compute 
their back-reactions on the geometry. The gravity waves associated with the collisions 
evolve with time and from the induced boundary metric, one can compute the energy 
momentum tensor of the field theory. Analyzing the time dependence of this stress 
tensor, one can learn about the evolution of the medium and subsequently account for 
the effects it has on the quark dynamics. 

Note that is different from Uk used in (l8|). 
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We have computed the fluxes and the form of the warp factor for our static dual 
geometry, but did not give exphcit expressions for the deformation of the internal flve 
dimensional metric to all order in QsNf. However, we have explicitly shown the Einstein 
equations that determine the form of the internal metric and using our ansatz in [1], 
one can in principle compute coefficients in the expansion of the internal metric to all 
orders in QgNf. Even without a precise knowledge of the internal geometry, we were 
able to extract crucial information about the dual gauge theory and formulated how the 
higher order corrections may enter into our analysis. Most of the calculation only relied 
on the warp factor and with the knowledge of its precise form, we were able to calculate 
thermal mass, drag and diffusion coefficients, rj/s and flnally free energy of QQ pair. 
For completeness of the supergravity analysis, we hope to compute the exact solution 
for the internal metric using the ansatz of |4] in our future work. 

In our computation of the heavy quark potential, we classifled the most general 
dual gravity that allows linear conflnement of quarks at large separation and small 
temperatures. We showed that if a gauge theory has dual gravity description and its 
effective degrees of freedom grows monotonically in the UV, it always shows linear 
confinement at large distances, as long as the dual warp factor satisfies a very simple 
relation given by (!52|) . Thus ( 152|) can be regarded as a sufficient condition for linear 
confinement of gauge theories with dual gravity. It would be interesting to study what 
are the general brane configurations that allow warp factors which satisfy ( l52l) and thus 
give rise to confining gauge theories. We leave it as a future direction to be explored. 

As the dual geometry incorporates features of Seiberg duality cascade, our 
construction is ideal for studying phase transitions. The gauge theories we studied 
have description in terms of gauge groups of lower and lower ranks. From the gravity 
dual analysis, by cutting the geometry at certain radial location and attaching another 
geometry up to infinity, we can construct gravity description for various phases of a 
gauge theory. Each phase will have different dual geometries attached in the large r 
region while the small r region will be common to all the theories. A flow from large r 
to small r geometry can be interpreted as 'flow' from an effective theory in the UV to 
another in the IR. From UV to IR, the different effective theories will describe different 
phases of the gauge theory and this can allow one to study the various phases of dense 
matter. Thus our construction is not only useful to analyze strongly coupled gauge 
theory, but also has potential for studying phase transitions in ultra dense medium and 
we hope to address this issue in the future. 

Acknowledgements 

Our work was supported in part by the Natural Sciences and Engineering Research 
Council of Canada, and in part by the Office of Nuclear Science of the US Department 
of Energy under grant No. DE-FG02-93ER40764. 



25 



References 



[1] K. Adcox et ai, (PHENIX), Nucl. Phys., A757, 184 (2005); J. Adams et at, (STAR), Nucl. Phys., 
A757, 102 (2005); B. B. Back et al, Nucl. Phys., A757, 28 (2005); I. Arsene et at, (BRAHMS), 
Nucl. Phys.,A757, 1 (2005). 

[2] See, for example, B. Schenke, larXiv: 1106.6012 [nucl-th]. and references therein. 

[3] M. Mia, K. Dasgupta, C. Gale and S. Jeon, Nucl. Phys. B 839, 187 (2010)t arXi?:0902.154Q [hcp-th]. 

[4] M. Mia, K. Dasgupta, C. Gale and S. Jeon, Phys. Rev. D 82, 026004 (2010) |arXiv:1004:0387l 
[hep-th]]; Phys. Lett. B 694, 460 (2011) [arXiv:1006.0055 [hep-th]]. 

[5] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 (1999)] 



arXiv:hep-th/9711200|; 



[6] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998) arXiv:hep-th/9802l50 ; S. S. Gubser, 
I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998) arXiv:hep-th/9802109 . 

[7] L R. Klebanov and M. J. Strassler, JHEP 0008, 052 (2000) |arXiv:hep-th/0007191 . C. Vafa, J. 
Math. Phys. 42, 2798 (2001) |arXiv:hep -th/0008142 . J. M. Maldacena and C. Nunez, Phys. Rev. 
Lett. 86, 588 (2001) ■arXrv:hep-th/0008001 . I. R. Klebanov and A. A. Tseytlin, Nucl. Phys. B 
578, 123 (2000) ,arXiv:hep-th/0002159, . L R. Klebanov and A. A. Tseytlin, Nucl. Phys. B 578, 
123 (2000) 

[8] P. Ouyang, Nucl. Phys. B 699, 207 (2004) arXiv:hep-th/0 311084] ; H. -Y. Chen, P. Ouyang, 

G. Shiu, JHEP 1001, 028 (2010). [arXiv:0807.2428 [hep-th]]. 
[9] C. Vafa, Nucl. Phys. B469, 403-418 (1996). hep-th/9602022 . 
[10] A. Karch and E. Katz, JHEP 0206, 043 (2002) [arXiv:hep-th/020523^ . 



[11] T. Sakai and S. Sugimoto, Prog. Theor. Phys. 113, 843 (2005) [arXiv:hep-th/0412141 . 

[12] O. Aharony, A. Buchel and A. Yarom, Phys. Rev. D 72, 066003 (2005) |arXiv:hep- th/0506002| ; 

JHEP 0611, 069 (2006) [arXiv:hep-th/060 8209 ; O. Aharony, A. Buchel and P. Kerner, Phys. 

Rev. D 76, 086005 (2007) |arXiv:0706.1768 [hep-th]]. 
[13] R. Casero, C. Nunez and A. Paredes, Phys. Rev. D 73, 086005 (2006) [arXiv:hep-th/060202r ; 

F. Bigazzi, A. L. Cotrone and A. Paredes, JHEP 0809, 048 (2008) arXiv:0807.0298 [hep-th]]; 

G. Bertoldi, F. Bigazzi, A. L. Cotrone and J. D. Edelstein, Phys. Rev. D 76, 065007 (2007) 
|a£Xiv:heH;h7070222f ; A. L. Cotrone, J. M. Pons and P. Talavera, JHEP 0711, 034 (2007) 
farXiv:0706.2766 [hep-th]]. 

[14] K. Dasgupta, G. Rajesh and S. Sethi, JHEP 9908, 023 (1999) 
K. Dasgupta, D. P. Jatkar and S. Mukhi, Nucl. Phys. B 523, 465 (1998) 



arXiv:hep-th/9908088 



arXiv:hep-th/9707224 



[15] S. B. Giddings, S. Kachru and J. Polchinski, Phys. Rev. D 66, 106006 (2002) 
arXiv:hep-th/0105097 . 



[16] M. J. Strassler, [arXiv:hep-th/0505153[ |arXiv:hep-th/0309"l49 



[17] K. Dasgupta and S. Mukhi, Nucl. Phys. B 551, 204 (1999) |arXiv:hep-th/9811139"; JHEP 9907, 
008 (1999) [arXiv:hep-th/9904131 ; K. Dasgupta, K. Oh and R. Tatar, Nucl. Phys. B 610, 331 
(2001) [arXiv:hep-th/0105066 ; JHEP 0208, 026 (2002) [arXiv:hep-th/0106040 ; K. Dasgupta, 
K. h. Oh, J. Park and R. Tatar, JHEP 0201, 031 (2002) [arXiv:hep-th/0110050 . 

[18] S. S. Gubser, C. P. Herzog, L R. Klebanov and A. A. Tseythn, JHEP 0105, 028 (2001) 
arXiv:hep-th/0102172 . A. Buchel, C. P. Herzog, I. R. Klebanov, L. A. Pando Zayas and 
A. A. Tseytlin, JHEP 0104, 033 (2001) |arXiv:hep-th/0102105). 

[19] M. Mia, arXiv:1012.2881 [hep-th]. 

[20] F. Chen, K. Dasgupta, P. Franche, M. Mia and S. Vaidya, To Appear. 



[21] M. Henningson and K. Skenderis, JHEP 9807, 023 (1998) |ar Xiv:hep-th/9806087| ; K. Skenderis, 
Class. Quant. Grav. 19, 5849 (2002) |arXiv:hep-th70209067] ; S. de Haro, S. N. Solodukhin and 
K. Skenderis, correspondence," Commun. Math. Phys. 217, 595 (2001) |arXiv:hep-th/0002230''; 
K. Skenderis, Int. J. Mod. Phys. A 16, 740 (2001) ^arXiv:hep-th/0010138 ; A. Karch, 



A. O'Bannon and K. Skenderis, JHEP 0604, 015 (2006) |arXiv:hep-th/0512125 



[22] G. W. Gibbons and S. W. Hawking, Phys. Rev. D 15, 2752 (1977). 



26 



[23] P. M. Chesler and L. G. Yaffe, Phys. Rev. Lett. 99, 152001 (2007) |arXiv:0706.0368l [hep-th]]; 

Phys. Rev. D 78, 045013 (2008) 'arXiv:0712,0050i [hep-th]] ; P. M. Chesler, K. Jensen, A. Karch 

and L. G. Yaffe, arXiv:0810.1985 [hep-th]. 
[24] Y. Kats and P. Petrov, JHEP 0901, 044 (20091. [arXrv:0712.0743l [hep-th]. 

[25] M. Brigante, H. Liu, R. C. Myers, S. Shenker and S. Yaida, Phys. Rev. D 77, 126006 (2008) 

;arXiv:0712.0805 [hep-th]]; 
[26] M. Brigante, H. Liu, R. C. Myers, S. Shenker and S. Yaida, Phys. Rev. Lett. 100, 191601 (2008) 

|arXiv:0802.33T8l [hep-th]]; 
[27] A. Dobado, F. J. Llanes-Estrada and J. M. T. Rincon, AIP Conf. Proc. 1031, 221 (2008) 

;arXiv:0804.2601 [hep-ph]]; 
[28] A. Buchel, R. C. Myers, M. F. Paulos and A. Sinha, Phys. Lett. B 669, 364 (2008) |arXiv:0808.T837l 

[hep-th]]; 

[29] I. P. Neupane and N. Dadhich. [arXiv:0808. 19191 [hep-th] : 

[30] X. H. Ge, Y. Matsuo, F. W. Shu, S. J. Sin and T. Tsukioka, JHEP 0810, 009 (2008) 

|arXiv:0808.2354l [hep-th]]; 
[31] A. Adams, A. Maloney, A. Sinha and S. E. Vazquez, JHEP 0903, 097 (2009). [arXiv:0812.0166l 

[hcp-th]. 

[32] A. Buchel, R. C. Myers and A. Sinha, JHEP 0903, 084 r2009'l. [arXrv:0812.2521l [hep-thj. 

[33] J. M. Maldacena, Phys. Rev. Lett. 80, 4859 (1998) |arXiv:hep -th79803002 . 

[34] J. P. Blaizot and E. lancu, Phys. Rept. 359, 355 (2002) 'arXiv:hep-ph/010lT03'. 

[35] J. Casalderrey-Solana and D. Teaney, JHEP 0704, 039 (2007) arXiv:hep-th/0701123 . 

[36] S. S. Gubser, NucL Phys. B 790, 175 (2008) |arXiv:hep-th/06 12 143j ! ' 

[37] B. R. Greene, A. D. Shapere, C. Vafa and S. T. Yau, NucL Phys. B 337, 1 (1990). 

[38] G. D. Moore and D. Teaney, Phys. Rev. C 71, 064904 (2005) arXiv:hep-ph/0412346l . 

[39] M. Mia and C. Gale, To Appear. 

[40] S. S. Gubser, Phys. Rev. D 74, 126005 (2006) |arXiv:hep-th/0605182| ; 

[41] C. P. Herzog, A. Karch, P. Kovtun, C. Kozcaz and L. G. Yaffe, JHEP 0607, 013 (2006) 

arXiv:h ep-th/0605158] . 
[42] K. G. Wilson, Phys. Rev. D 10, 2445 (1974). 
[43] L. Susskind, Phys. Rev. D 20, 2610 (1979). 
[44] A. M. Polyakov, Phys. Lett. B 72, 477 (1978). 
[45] C. Gale and J. 1. Kapusta, Phys. Lett. B 198 (1987) 89. 
[46] T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 

[47] S. Digal, O. Kaczmarek, F. Karsch and H. Satz, Eur. Phys. J. C 43, 71 (2005) 
;arXiv:hep-ph/0505193j ; 

[48] M. Doring, S. Ejiri, O. Kaczmarek, F. Karsch and E. Laermann, PoS LAT2005, 193 (2006) 
arXTvyhep-lat/0509150 . 

[49] O. Kaczmarek, F. Karsch, F. Zantow and P. Petreczky, Phys. Rev. D 70, 074505 (2004) [Erratum- 

ibid. D 72^059903 (2005)] arXiv:hep- lat/0406036 ; 
[50] K. Petrov, [arXiv:hep-lat /06l604lf . 

[51] Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz and K. K. Szabo, Nature 443, 675 (2006) 
;arXiv:hep-lat/06110l4 . 

[52] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. B 643, 46 (2006) 
arXiA^hep-lat/0609068 

[53] Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, S. D. Katz, S. Krieg and K. K. Szabo, JHEP 0906, 088 

(2009) arXiv:0903.4155 [hep-lat]]. 
[54] N. Brambilla, A. Pineda, J. Soto and A. Vairo, NucL Phys. B 566, 275 (2000) 

arXiv:heiriph/9907240|; Rev. Mod. Phys. 77, 1423 (2005) [arXiv:hep-ph/ 0410047| ; 
[55] N. Brambilla et al. [Quarkonium Working Group], arXiv:hep-ph/0412158; 

[56] N. Brambilla, J. Ghiglieri, A. Vairo and P. Petreczky, Phys. Rev. D 78, 014017 (2008) 
|arXiv:0804.0993l [hep-ph]] . 



27 



[57] J. Polchinski and M. J. Strassler, Phys. Rev. Lett. 88, 031601 (2002) [arXiv:hep-th/0109174 

JHEP 0305, 012 (2003) |arXiv:hep-th/0209211 . 
[58] H. Boschi-Filho and N. R. F. Braga, JHEP 0305, 009 (2003) |arXiv:hep-th/0212207^ ; Eur. Phys. 

J. C 32, 529 (2004) [arXiv:hep- th/0209080| ; 
[59] H. Boschi-Filho, N. R. F^ Braga and C. N. Ferreira, Phys. Rev. D 73, 106006 (2006) 

[Erratum-ibid. D 74, 089903 (2006)] l arXiv:hep-th/0512295j ; Phys. Rev. D 74, 086001 (2006) 
arXiv:hep-th/0607038. 

[60] C. A. Bahon Bayona, h7 Boschi-Filho, N. R. F. Braga and L. A. Pando Zayas, Phys. Rev. D 77, 

046002 (2008) [arXiv:0705.1529 [hep-th]] 
[61] M. Panero, Phys. Rev. Lett. 103, 232001 (2009) arXiv:090L^T9" [hep-lat]] . 



[62] S. J. Rey S. Theisen and J. T. Yee, Nucl. Phys. B 527, 171 (1998) arXiv:hep-th/9803135]; 

[63] S. J. Rey and J. T. Yee, Eur. Phys. J. C 22, 379 (2001) |arXiv:hep-th/9803001 

[64] A. Brandhuber, N. Itzhaki, J. Sonnenschein and S. Yankielowicz, Phys. Lett. B 434, 36 (1998) 

arXiv:hep-th/9803137 ; from supergravity" JHEP 9806, 001 (1998) arXiv:hep -th /9803263 
[65] D. J. Gross and H. Ooguri, Phys. Rev. D 58, 106002 (1998) [arXiv:hep-th/9805129j . ' 
[66] F. Bigazzi, A. L. Cotrone, A. Paredes and A. Ramallo, Fortsch. Phys. 57, 514 (2009) 

arXiv:0810,.5220 [hep-th]]; JHEP 0903, 153 (2009) [arXivi0812,3399 [hep-th]]. 
[67] E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane and T. M. Yan, Phys. Rev. Lett. 36, 500 

(1976); Phys. Rev. D 17, 3090 (1978) [Erratum-ibid. D 21, 313 (1980)]; Phys. Rev. D 21, 203 

(1980). 

[68] O. Andreev and V. I. Zakharov, Phys. Rev. D 74, 025023 (2006) |arXiv:hep-ph/0604204] ; 
Phys. Lett. B 645, 437 (2007) |arXiv:hep-ph/0607026|; JHEP 0704, 100 (2007) 



arXiv:hep-ph/0611304j ; Phys. Rev. D 76, 047705 (2007) |arXiv:hep-ph/0703010| . 



